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We consider the dynamics of a disordered p-spin model with p — 2, analyzing the dynamics within 
Ruelle's thermodynamic formalism, We use an indicator of the dynamical activity to construct the 
relevant dynamical Gibbs ensemble. We prove that the dynamics in the low-temperature (spin 
glass) phase of the model take place at a second-order phase transition between dynamically active 
and inactive trajectories. We also show that the same behaviour is found in a related model of a 
three-dimensional ferromagnet. 
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I. INTRODUCTION 

Glassy systems are characterised by their dynamical 
properties: at their glass transitions, they fall out of equi- 
librium on experimental time-scales, and exhibit aging 
phenomena. As the glass transition is approached, their 
relaxation times increase in a super- Arrhenius fashion 
and the decay of their equilibrium dynamical correlation 
functions is slower than exponential [l| . In the last fifteen 
years, several approaches have been developed to recover 
theoretically the experimental, out-of-equilibrium results 
(see, for example, [2|)- 

One of the most striking experimental features of a 
glass-forming liquid is that the increase in relaxation time 
near the glass transition does not seem to be accompanied 
by any significant changes in the liquid structure. How- 
ever, experiments and computer simulations [1, 0, Q 
both indicate that a dynamical length scale is growing 
as the glass transition is approached. That is, glassy 
materials are made up of active and inactive regions 
of space-time, namely dynamical heterogeneities. Based 
upon these results, the idea that the glassy properties of 
a system arise directly from their dynamical heterogene- 
ity was developed in 0, [1] , stimulating further theoretical 
understandings of glassy dynamics. 

Spin glasses are magnetic spin systems which exhibit 
several features in common with glass-forming liquids. 
They are modelled by spins with quenched random inter- 
actions between them, and have been extensively inves- 
tigated both experimentally and theoretically (see for 
a review). The purpose of this paper is to show that the 
glassy dynamics of a particular spin glass model can be 
understood in terms of the histories it follows in configu- 
ration space. To this end, we employ the thermodynamic 
formalism of histories, developed by Ruelle and coworkers 
within the framework of dynamical systems theory, 
and summarized in [TTj in the context of Markov dynam- 
ics. While equilibrium statistical mechanics is concerned 
with fluctuations in the configuration space of the sys- 
tem, Ruelle's formalism focuses on the trajectories (his- 
tories) by which the system evolves through configuration 
space. The method has been applied recently to kinet- 



ically constrained models of glass-formers [iJI and to a 
Lennard Jones binary mixture [isj . Both these studies 
distinguish active and inactive histories of the systems, 
according to the range of configuration space visited dur- 
ing the history. In the kinctically constrained models, it 
was proven that the active and inactive histories form dis- 
tinct populations. In the language of the thermodynamic 
formalism, they are separated by a first-order phase tran- 
sition in trajectory space. In Refs. [H, [l^j it was argued 
that the heterogeneous dynamics of those models is in- 
trinsically linked to this transition. 

In this article, we focus on a soft p-spin model with 
p = 2, whose static and dynamic properties can be stud- 
ied analytically. The p oo limit of the p-spin model, 
namely the disordered Random Energy Model, was re- 
cently shown to possess a connection between the ac- 
tivity of the histories it follows and the dynamical het- 
erogeneities in its glassy phase [ij]. (Of course, in p- 



spin models with infinite-ranged interactions, the dy- 
namical correlation lengths associated with dynamical 
heterogeneity are ill-defined. However, the presence of 
large dynamical fluctuations in these mean-field models 
is naturally linked to dynamical heterogeneity in their 
finite-dimensional counterparts. We show that the p = 2 
spin-glass model is closely related to a three-dimensional 
ferromagnet in which such length scales can be calcu- 
lated.) In both the spin glass and ferromagnetic models, 
we demonstrate the existence of a phase transition in tra- 
jectory space, offering further evidence that these phase 
transitions are very generally associated with glassy sys- 
tems. 

The outline of the paper is as follows: in Section |TT] 
we describe the models we will consider and the meth- 
ods that we will use. In section IIIII we construct a 'dy- 
namic phase diagram' that describes the behaviour of the 
system in trajectory space. We interpret our results in 
Sec. IIVI discussing the links between the large deviations 
that we have derived and the more familiar features of 
the soft spin models, and identifying directions for fur- 
ther study. 
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II. MODEL AND FORMALISM 

A. Spin glass model 

We consider a system of N continuous spins ai whose 
Hamiltonian Ti is given by: 



2N 



u 
N 



(1) 



Here, (3 = 1/T as usual, where T is the temperature and 
we have set Bohzmann's constant to unity, we take m > 0, 
and the random couplings arc Gaussian distributed 



/27r 



exp 




(2) 



The role of the term proportional to u > is to sup- 
press configurations with extreme values of the spins. 
The model is similar to the p-spin models discussed 
in [l^, [31 ■ Like the spherical p-spin model of [l^ and 
in contrast to that of [l^, the model under considera- 
tion here can be solved exactly. However, we use the u- 
tcrm instead of a spherical constraint since it facilitates 
studies of large-deviations of the activity. (In particu- 
lar, we note that fluctuations of extensive quantities re- 
quire careful treatment in spherical models [17| , in which 
fluctuations of the Lagrange multiplier for the constraint 
must be considered.) In any case, all of these soft p-spin 
models exhibit flnite-temperature 'glass transitions' at 
which ergodicity is broken [T^, [H, [l^, [131 ■ Connections 
with the structural glass problem have been discussed 
in [H, [m, [Hi. The case of p = 2 differs from that of 
p > 3 in that correlation functions can be obtained ex- 
actly from the properties of large random matrices [ist . 

In this article, we will employ the functional-integral 
formalism of [23j . To this end, we endow the spin system 
with a relaxational dynamics, 



dtCTi = - 



6f3H 
S<J,{t) 



(3) 



where the t^^'s are independent white Gaussian noises 
with variance 2. Following [3 [13] ; will prove use- 
ful to resort to the basis which diagonalizes the matrix of 
exchange couplings. The eigenvalues { Jp}/j=i....,jv of the 
N X N matrix {Jij)ij=i^...^N are distributed according to 
Wigner semi-circle law. 



(4) 



Denoting by 0^ the spin coordinates in the basis that in 
which (Jij) is diagonal, the Hamiltonian simpliflcs into 



(5) 



and the equation of motion now reads 



where the 77^ 's are independent white Gaussian noises 
with variance 2. 



B. Ferromagnetic model 

It been remarked that the p = 2 spin glass resembles a 
ferromagnet 'in disguise' [13]. To illustrate this, we also 
consider a ferromagnetic model whose Hamiltonian TYfm 
is given by 



FM 



(7) 



where the flrst sum runs over nearest neighbours on a d- 
dimensional (hyper)-cubic lattice, but the u-term retains 
interactions between all sites. (Thus, the model contains 
infinite-ranged couplings, as in the spherical ferromag- 
net.) The analogues of the co-ordinates (j)^ in this model 
are the Fourier transformed spin co-ordinates (jjk where 
k = (fci, . . . , kd) is the wave vector. The eigenvalues of 
the matrix coupling the spins are Ek = X]r=i ^r- The 
resulting equation of motion is then 



d_ 
dt 



^Ei'^'= 



Vkit) (8) 



where the rjk arc independent Gaussian noises as before. 

In d = 3 the distribution of the eigenvalues Ek is 
p{Ek) ~ (27r^)^^ when |fc| is small. Similarly, 
in the spin glass of ^ , the density of eigenvalues scales 
as p{J^) ~ TT^^y^{2 — Jf^) for close to 2. We will find 
that the phase transitions in the models depend on the 
scaling of the eigenvalue density near these points, and 
hence that phase transitions in the d = 3 ferromagnet 
and the p = 2 spin glass are related to each other, and 
have the same scaling exponents. 



C. Symmetry-breaking fields 

Below its transition temperature, the ferromagnetic 
model spontaneously breaks the global symmetry ai 
—<Ti. To clarify the behaviour in the ordered phase, 
it is convenient to introduce a magnetic field: we take 
''^FM ^ "^FM ^ '^X^i "'i ™ ©• ^^'^ equation of motion 
becomes 



d 1 

—6k = l3{Ek(f>k + hSkfiVN) + 4u0fc— V \(t>k\'^ + T]k{t) 
at Jy ^-^ 

k 

(9) 

In the presence of this field, the magnetisation A4{t) = 
N~^^-ai = N~^/'^(l)k=o acquires a finite expectation 
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value: in the ordered phase, the magnetisation remains 
finite even in the Umit of small h, with a first-order phase 
transition at h = 0. 

In the spin glass model, the low temperature phase 
breaks the same global symmetry, but the order parame- 
ter is not the magnetisation. Instead, the co-ordinate t/);^ 
corresponding to the largest eigenvalue of Jij becomes 
macroscopically occupied. We assign the label = 
to this eigenvector. Then, by analogy with the ferromag- 
netic case, we can introduce an analogous 'staggered field' 
to the model of (H]). The equation of motion becomes 



d 



1 V-^' 

N ^ 



(10) 

The analogue of the magnetisation M (t) is the staggered 
magnetisation A4s{t) = N~^/'^(j)^^Q. In the low temper- 
ature phase of the p = 2 spin glass, the expectation of 
A^s(^) tends to a finite value as hg tends to zero, with a 
first-order phase transition at = 0. 



D. Thermodynamic formalism 

Ruelle's thermodynamic formalism involves a statisti- 
cal mechanical analysis of the trajectories that a system 
follows through configuration space. Let a history be 
a particular time-realization that the system has visited 
over a given time interval. Consider an ensemble of histo- 
ries constructed by fixing their dynamical activity K{t). 
Here, the dynamical activity is a history-dependent ob- 
servable, extensive both in space and time, expressing 
the amount of activity within the history. In a typical 
inactive history, the spins remain frozen in a locally or- 
dered state; in an active history, spins fluctuate randomly 
between up and down states. For both spin-glass and fer- 
romagnetic models, a simple local observable consistent 
with this definition of activity is 

^^'W = -^E Tdt'^lw^-^E Tdi^^w (11) 

^ ^ ^ ^ JO 

Histories with K close to are the most active ones and 
are typically associated with disordered states; histories 
with large negative K are inactive and are associated 
either with local or global ordering. 

While an ensemble of trajectories with fixed K is nat- 
ural from a physical point of view, our theoretical meth- 
ods require a change of ensemble, to one in which the 
average activity is fixed. (To draw an analogy with equi- 
librium statistical mechanics, we are transforming from 
a microcanonical to a canonical ensemble.) To fix the 
average activity, we apply a field s that is conjugate to 
K{t). While we are as yet unable to endow s with an 
experimentally-realizable physical meaning, ensembles of 
histories with finite s provide a valuable theoretical tool, 
which allow us to probe the histories that the system fol- 
lows. The ensemble with s = is simply the (unbiased) 



ensemble of trajectories for the system: ensembles with 
s > are less active than the unbiased ensemble while 
those with s < arc more active. 

In the following, we will evaluate the partition function 



Z{s,t) - 



(12) 



which is simply the generating function for the activity. 
Here and throughout, we use (•)o to denote an average 
of the (unbiased) rclaxational dynamics over all possible 
time realizations, which means an average over the noises 
rji . We also consider averages of a generic observable A in 
the biased ensemble parameterised by s, which we write 
as {A)s = limt_oo^~'(s,0(^e-^^)o 
We also define a dynamical free energy 



V'(s) 



lim 

t — »oo 



\nZ{a,t) 



It follows that 



(AO. 



^Nt 



d7 



(13) 



(14) 



With these definitions, singularities in tl){s) are the dy- 
namical phase transitions of the system. Discontinuities 
in the derivatives of ^^{s) will correspond to phase tran- 
sitions between active and inactive phases. By analogy 
with equilibrium statistical mechanics, transitions with a 
jump in {K) s are termed 'first-order' or 'discontinuous'; 
otherwise the transition is termed 'continuous'. More 
specifically, if {K)s is continuous and the second deriva- 
tive is discontinuous then we refer to the transition as 
'second-order'. 



E. Functional integral formulation 

To evaluate dynamical observables such as (i^)s, we 
use the Janssen-De Dominicis functional-integral formu- 
lation [2^, as in earlier studies such as Ref. [l^. The 
rclaxational dynamics for the spin 4>fj,{t) is given by 
Using a functional integral representation the dynam- 
ical, s-dependent partition function introduced in (TH 
becomes 



Z{s,t) = / V(j>V(l)exp 



- / di'L(i') 



(15) 



where, omitting time-dependence for brevity. 



L = 



4u ' 



N ^ 



(16) 

(We consider the spin-glass model of ((T|) and we have set 
the staggered field = 0; other cases will be discussed 
below.) 

Due to the infinite-ranged interactions in the term pro- 
portional to u, this model may be reduced to a quadratic 



4 



form. Details arc given in appendix[Xl The result is that 
the partition function becomes 



Z{s,t) 
with 



Vcj/Dcj) exp 



^uxxNt- I dt'iaux(i') 




(17) 



^aux = ^ 0/.(9t' + 4UX - PJfi 



(f-4.x) 



4wx ) <f>l 
(18) 

where the parameters x X must be determined self- 
consistently, through 



X= dJ^p(J^) 



1 



v/(4ux-/3J^)2-2(s-8ux) 



(19) 



and 
2X 



1= / dJ^p(J^)- 



4ux - PJ^, 



V(4«x - /3 J^)' - 2(s - 8ux) 

(20) 

In addition, it follows from the definition of x that 



X = -2(/^). 



(21) 



so that solving the self-consistency equations ([19]) and 
([20]) leads directly to the activities of the relevant phases. 
Finally, we note that the derivative 



dx 
ds 



(22) 



gives the fluctuations of the activity. 

Physically, we have shown that the dynamical correla- 
tion functions of the original model (H]) are the same as 
those of the auxiliary quadratic system (|17p . Noting that 
p{J^) is finite only for —2 < < 2, the integrals in ^ 
and (Uni) are well-defined as long as 



(4mx - 2/3)2 - 2(s - 8ux) > 



(23) 



As long as this condition is fulfilled, then the system is 
in a paramagnetic disordered phase. On the other hand, 
if the denominator of vanishes at = 2 then the 
mode associated with this eigenvalue may become macro- 
scopically populated. 



F. Symmetry-breaking fields 
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FIG. 1: 'Phase diagram' associated with the dynamic free en- 
ergy 'i/'(s) of the spin glass model. The critical point of the 
model is at s = and T = Tc. The heavy solid line is a sec- 
ond order phase boundary between ordered and disordered 
phases. The dashed line is a crossover within the disordered 
phase. In the high temperature (normal) regime then the re- 
sponse to a staggered field hs is positive; the low temperature 
(anomalous) regime is characterised by a negative response to 
this field. 



respectively. Evaluating these expectation values in the 
auxiliary model, we arrive at 



Ills ~ ph 
iris = Ph 



Aux - 2/3 



(4ux - 2/3)2 - 2{s - 8ux) 

s - Sux 

(4ux - 2/3)2 - 2(s - 8ux) ' 



(24) 
(25) 



Self-consistency in the presence of the field leads to mod- 
ified saddle point equations for x and x- 



X = + 
2x + 1 = 2msfrL. 



(26) 



dJ^p{J^){AuX'' PJf,) 

^{Aux^pJ^Y-2{s~8ux)' 

(27) 



Finally we note that while we have considered the spin- 
glass model of ([T]) throughout Sees. Ill El and III F( the 
equations for the ferromagnetic model can be obtained 
by applying the simple replacement (/x, J^, mg, rn-s) 
{k,Ek,m,m) throughout these sections, where m = 

iV-l/24=0- 



As discussed in Sec. Ill C\ it is also useful to consider the 
effects of symmetry-breaking fields hg and h on these sys- 
tems. Following the analysis of the previous section, the 
symmetry-breaking fields lead to linear terms in iaux- In 
general, the fluctuating magnetisation Ms = N^^/'^(t>^=o 
and its response fleld Ms = N~^/^(j}f^^o both have fi- 
nite expectation values which we denote by mg and fhs 



III. DESCRIPTION OF THE PHASE DIAGRAM 

A. Overview 

The dynamical phase diagram for the p = 2 spin glass 
of IH) is shown in Fig. [1] for /ig = 0. The phase diagram 
of the ferromagnetic model of ([7]) at ft, = has the same 
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form. The axis s = corresponds to the unbiased re- 
laxational dynamics of ([31) ■ On this axis, the spin-glass 
system has a second-order phase transition at 



1 



2u 



(28) 



Below Tc, the mode with lowest eigenvalue is macroscop- 
ically populated, and the system is ordered: the (stag- 
gered) magnetisation remains finite as the (staggered) 
field is reduced to zero. 

We will show that the effect of positive s is to pro- 
mote ordering in the system, consistent with the expec- 
tation that ordered phases are less active that disordered 
ones. As s increased from zero, the second-order transi- 
tion between active and inactive phases moves to a higher 
temperature: we have an increasing function Tc{s), with 
7c(0) = Tc being the thermodynamic transition temper- 
ature at s = 0. 

The effect of negative s is to reduce ordering in the sys- 
tem, thus increasing the activity: at temperatures above 
Tc, the dynamical free energy ipis) has no singularities 
for s < and K increases smoothly as s is decreased from 
zero. We also find that no ordered phases are possible for 
s < 0: the condition ((23)) is always satisfied when and nig 
vanishes in the limit of small hs- Thus, for < T < Tc, 
behaviour of the model as s — > 0^ coincides with the or- 
dered phase that is found at s = 0, but the behaviour 
is different for all s < 0. This signals the presence of a 
phase boundary at s = 0. The same effect is observed in 
the ferromagnetic model of ([7]). 

We now show how this phase diagram is obtained 
from the solutions to the self-consistency equations (fTO)) 
and (PO)) . We calculate the saddle point average x as a 
function of s and the other parameters of the model: this 
gives the activity of the phases of the model through (PT|) . 
In addition, we also calculate the (s-dependent) stag- 
gered magnetisation mg which gives additional insight 
into the phases of interest. 



B. Unbiased dynamics (s = 0) 

We begin with the unbiased (s — 0) behavior of the 
model, the derivation of which is identical to that of 
Kosterlitz et al [T§| for the spherical version of this model. 
Equations (|I9p and (PO)) can be solved at s = with the 
result that x — ^i as required by causality, and that 



AT ^ 



I 



(29) 



If X > ^ then we can approximate the sum over the 
eigenvalues by an integral over the distribution p{J^), 
with the result 



^ (2ux - V4«2x2 - 



(30) 



the solution of which is given by 
1 



X ■ 



(31) 



where /?c = 1/Tc = ^/2u consistent with 

However, for (3 > Pc, the mode with = 2 becomes 
macroscopically occupied, as described above. We there- 
fore have 



X 



2u' 



(32) 



While the integral in (fTQ]) is formally undefined, the or- 
dered phase can be studied either by introducing a finite 
staggered field as discussed in Sec. Ill Fl or by solv- 
ing (fT9|) to 0{1/N). In either case, the staggered mag- 
netisation m„ at zero field is 



Too 



T2 - r2 



T 



(33) 



We also note a property of the ordered phase that is 
peculiar to exactly soluble soft spins models such as the 
spherical model ^M,- The susceptibility associated with 
the magnetisation. 



to(2) = {{M{t)-m,f), 
— > for the unbiased dynamics (s 



(34) 



diverges as h. 

and all T < Tc- (This can be verified by evaluating ((/)q) 
in the auxiliary model.) This is in contrast to the usual 
situation in critical phenomena where to^^) jg finite for 
zero field and T < Tc, diverging only at the critical point. 



C. Ordered phase, s > 

We now turn to the ordered phase for positive s. To 
simplify the analysis, we introduce reduced variables: 



Xis) 
Yis) 



AuTx{s) 
2T^\s- 



Suxis)] 



(35) 
(36) 



where we explicitly indicate the s-dependence of x and 
X- We see that X{s) = A*: physically, we identify the 
quantity (— A*(/)^) in (jA3| as the constraint force on mode 
fi that arises from the w-term in ^ and suppresses con- 
figurations with extreme values of 0^. We also identify 
Y{s) as a renormalised field s for the auxiliary system. 
It is easily verified that while the dynamical free energy 
■0 depends on four parameters (s, /3, u, h^), the properties 
of the auxiliary model depends only on (Y, (3, X, hg). Our 
strategy will be to determine properties of the auxiliary 
model in terms of X and Y and then to find the relations 
between {X, Y) and the bare parameters of the model. 

In the presence of a staggered field hg, we have 
from dMl) that 



X- 



Y 



(37) 
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For y > 0, spontaneous symmetry breaking occurs if the 
denominator vanishes as ft., — *■ 0, as 



X 



Y + 0{K) 



(38) 



Working at smaU /ig, we then take the zeroth order terms 
in (|26p . arriving at 



Y 



AuT 



ml+T j dJp{J) 



1 



2(2- J)VF+ (2- J)2 



1/2 



(39) 

This allows us to obtain = ml + ^^Y^^"^ + 0{Y^/^). 

Finally, we must relate the renormalised field Y to the 
bare field s. Again working at zeroth order in hg, ([27]) 
becomes: 



2s - y/?^ 

8u 



PmiVY - 1+ 



dJp(J) 



2-J + VY 



2(2- J)VF+(2-J)2 



1/2 



(40) 



Taking s small and positive, the solution has small pos- 
itive Y. More specifically, the first term on the right 
hand side of ((40l) dominates as y 0, leading to 
VY = 2(ff^-2M) • [The same result can be obtained by 
working at /ig = and considering carefully the limit of 
large-iV. The analogue of dSS]) is X = 2 + Vy + 0{1/N) 
and we allow for a finite values of rus and ifis when solv- 
ing (pS)) and (P7)) . The remainder of the analysis follows.] 
Taking everything together, in the limit hg —>■ 0+ and 
for < s ^ (T — Tc), the leading behaviour is 



s(3 



2u 8u(/32-2m) 



7Tf3y (3^~ 2u 



(41) 
(42) 



Physically, we can see that the s = axis in the phase 
diagram of Fig. [l]is singular, but that both x and d^/ds 
are finite, so that fluctuations of the activity K(t) remain 
finite as /i ^ 0. However, as for the case s = 0, the 
fluctuation wS"^^ diverges as /i ^ 0, for all cases where 
the spin-reversal symmetry is spontaneously broken. 



D. High temperature regime 

We now turn to temperatures above the critical tem- 
perature, /3 < /3c. We treat s perturbatively in and 
(PO)) . arriving at 



V4u - /32 8(2w-/32)^4u-/32 
noting also that 

2U-/32 



y = 2s 



/?2(4u-/32) 



(43) 



(44) 



We also notice that these solutions satisfy (X — 2)2 — y > 
at small s. 

As the temperature is lowered towards Tc, we see that 
dx/ds diverges. This again signals the second-order tran- 
sition to the ordered phase. Indeed, it can be shown that 
this transition to the ordered phase moves to a higher 
temperature for s > 0. At the critical point, we have 
{X — 2)2 = y and mg = rhg = 0. With these conditions, 
we can use (|39p and (|40p to derive the phase boundary 
for positive s. In the limit /3 /?~, the system is ordered 
for s > Sc, where 



Sc(/3) 



2^ 

3/3c 



(/3c - /3)= 



(45) 



which holds to leading order in /3c — /?■ This function 
gives the phase boundary in Fig. [1] and its inverse gives 
the function Tc(s) discussed above. 

This completes our analysis of the phase diagram for 
s > 0. There is an ordered phase separated from a para- 
magnet by a second-order phase transition. Loosely, the 
effect of positive s is simply to stabilise the ordered phase, 
so that spontaneous symmetry breaking takes place at a 
higher temperature. 



E. Anomalous paramagnetic regimes 

We now take s < but we remain in the low tem- 
perature regime with (3 > (ic- Working in terms of the 
reduced variables X,Y , we take y < so that we have 



(a: - 2)2 - y > 



(46) 



and the integrand of ([^D|) is finite for — 2 < < 2. To 
make progress with the integrals in (|19l20p : we define 



"^{X-J,)^~Y 



(47) 
(48) 



We will consider the limit Y , in which the solution 
to (HOD is a: ^ 2". The relevant limit is < {-Y) < 
(2 - X)2 < 1. Writing y' = -Y and x' = 2 - AT, and 
after some manipulations presented in appendix [Bl we 
obtain: 



V 



h = TT- Vx' In {Ax'-'/y') + 1 + + 0(x'i/2) (49) 



!(^)!^! + o(yV"i/2) + o(^'5/2) 



Btt 



The self-consistent equation (|19p takes the form 

^(2-a;') = l + — In— +0(1 , 
4m tt y 



(50) 



(51) 
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where we define o(l) terms as quantities which vanish for 



2/ < x'^ < 1 . RecaUing that X = 2 
have to leading order 



A, 2 

4w V iu^ln {^1/Y)J 



'iuTx, we 



(52) 



Then, substituting for x' in ([50]) and recalhng ((20|) . we 
have 



12^ W (-1/r) 



(53) 



Noting that y' = —2T^{s — 8ux) we see that when s — > 0, 
X » y', so that s ~ 8ux + o(l). This simphfies the 
expressions for x a-nd x which become: 



X 



2w 

s 

8w 



16^ 



(54) 



(55) 



which hold at leading order in s < and for (3 > (3c- 

We refer to the phase with s < and T < Tc as an 
anomalous disordered phase. To understand this ter- 
minology, it is useful to consider the linear response 
to the field /is, which is given by P^ . Since we have 
< -y < (2 - Xf with X >2, this reduces to 



2- X 



(56) 



We can see that the response to the staggered field is a 
staggered magnetisation in the opposite direction (a dia- 
magnetic response). In addition, it follows from (j54|l that 
drng/dft-s '-^ —\s\~'^/'^ for s 0^. That is, the diamag- 
netic response diverges. Clearly such a response would 
be impossible in the unbiased (s = 0) ensemble due to 
thermodynamic convexity arguments, but when consid- 
ering ensembles with finite s then such arguments do not 
apply. 

Finally, we consider the nature of the phase transi- 
tion between ordered and anomalous disordered phases. 
Comparing ([1T|) and ([51)) . we note that x is continuous 
at s = 0, and hence that the activity {K)s is continuous 
also. Thus, we identify a continuous phase transition at 
s = 0, consistent with Fig. [TJ Often, at continuous phase 
transitions, one may identify a path between the phases 
along which the free energy is analytic and which remains 
always near the critical point. (For example, in a ferro- 
magnet one can move between ordered and disordered 
phases by applying a small field h, decreasing the tem- 
perature and then removing the field.) However, in the 
transition considered here, the (staggered) magnetisation 
TOg is zero for s < but has a finite limit as s ^ 0^ . This 
seems to preclude such a route around the critical point. 
Further, evaluating (dx/ds) indicates that the fluctua- 
tions of the activity diverge as s ^ 0~ but remain finite 
as s ^ 0. This also indicates the absence of a path be- 
tween the phases that is continuous near the transition. 



We turn to this issue in the next section, where we also 
discuss the possibility of diverging length scales near this 
transition. 



F. Effect of the staggered field hs on disordered 
phases 

To understand the behaviour near this phase transition 
in more detail, we introduce a finite staggered field /is > 
and note that the crossover between normal (rus > 0) and 
anomalous (ms < 0) behaviour takes place at x = /3/2w. 
If we insist that x take this value, the self-consistency 
equation (QH) becomes 



1 



(57) 



For (3 < Pc (high temperatures), this equation has a 
unique solution for y < 0, which signals a crossover from 
normal to anomalous behaviour, at a value of Y that is 
independent of the field hs- On the other hand, in the 
low temperature regime (3 > Pc, (|57p has no solutions. 

One may verify that for large positive s, the solution 
of (pS)) has X ^ 0, while for small positive s we have 
from ([41]) that X > 2. Since ((57)) establishes that, for low 
temperatures, there there are no values of Y for which 
X ^ 2, it follows that the s-dependent value of X has a 
jump from a value greater than 2 to a value smaller than 
2. This is a discontinuous phase transition from normal 
to anomalous states. (The field hs is finite, so the concept 
of a spontaneously ordered state is not useful. However, 
the sign of the staggered magnetisation is positive in the 
normal state and negative in the anomalous one.) These 
arguments lead us to propose the qualitative phase dia- 
gram shown in Fig. [21 The presence of first-order tran- 
sitions at finite /is explains the unusual features of the 
order-disorder transition at s = /is = 0: there are indeed 
no continuous paths between the ordered and disordered 
states due to the first order transition at finite /is. 

The nature of the scaling behaviour near Tc in this 
model is clearly complicated, depending qualitatively on 
the order in which s, hs and T — Tc are taken to zero, and 
also on the signs of s and T — Tc- A detailed investigation 
of these finite-/is transitions is beyond the scope of this 
paper. However, we can conclude that the order-disorder 
transition at /is = s = is second-order in that xi^) is 
continuous, but that the spontaneous staggered magneti- 
sation TTig goes discontinously to zero at this transition. 
We again emphasise that all of this phenomenology is 
also present in the ferromagnetic model of ([7]), at least 
for d — 3. We expect the qualitative features to also be 
present in higher dimensions. 

In the ferromagnetic model, we can also consider the 
correlation lengths of the various phases. These appear 
through the fc-dependence of the fluctuations: 



T 



Y 



(58) 
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FIG. 2: (a) Proposed phase diagram for > 0. The sohd line is a first-order phase transition at which the staggered 
magnetisation and the activity X are discontinuous. It ends at a critical point at T = Tc and s — 0, but we note that ha is finite 
at this critical point. The dashed line is a crossover at which the linear response to a staggered field hs vanishes. The dashed 
line is independent of hs while the solid line approaches the s = axis as hs — > 0. (b) Behaviour as a function of the field hs in 
the high temperature regime. For negative s, there is a crossover from normal to anomalous response: at the crossover X = 2 
and Y satisfies (|57p . For positive s, there is a critical point at Sc, with spontaneous symmetry breaking for s > Sc. Close to 
Tc, the critical value of Sc is given by (|45p . (c) Behaviour as a function of hs in the low temperature regime. The solid line for 
s < is a first-order phase transition between states with positive and negative response to the staggered field, while the line 
for s > is the usual first-order transition between spontaneously-ordered states. 



Recall that < d with equality if fc = 0. 

Two cases are of interest. Firstly, if the denominator 
vanishes at fe = as in the ordered states, the fluctu- 
ations of the spontaneous magnetisation diverge, as de- 
scribed above. Secondly, if X < d as in the anomalous 
phase then S{k) has a peak at a finite wave- vector fe* for 
which Ek<- ~ X. Wc interpret l/|fc*| as a characteristic 
length scale for structures within this phase. It is inter- 
esting to note that this length scale diverges as s ^ 0" 
in the anomalous phase, and that this is accompanied by 
a divergence in the fluctuations of the spontaneous mag- 
netisation [it may be easily verified that 5'(0, s — *■ 0~) is 
divergent since X d~ and y — > 0~ in this limit, by 
analogy with the p = 2 spin glass] . 



IV. INTERPRETATION 

We have considered in some detail the large deviations 
of the dynamical activity in two soft-spin models. We end 
with a comparison with previous studies and with some 
comments on the relation between the large deviations 
that we studied and the phase behaviour of the models. 

Based on the close relationship between the ferromag- 
netic and spin glass models, the form of the phase dia- 
gram in Fig. [1] is perha ps n ot too surprising: a similar re- 
sult was found in Ref. Tl'| for the infinite-ranged (mean- 
field) Ising model. The equilibrium critical point leads to 
a dynamical phase transition at s = below the critical 
temperature. However, instead of being first order as in 
the fully-connected model, the phase transition we have 
found is second order. Unlike the fully-connected mod- 
els, the ferromagnetic model of ([7]) has several diverging 
length scales, although the presence of a diverging length 
scale throughout the ordered phase may be a peculiarity 
of our particular model. In any case, the diverging length 



scale as s — > 0~ within the anomalous phase seems to be 
a new feature that merits further investigation. 

In particular, the existence of the anomalous phase 
seems to be linked to the existence of aging/coarsening 
solutions to the relaxational dynamics of these models. 
These solutions are characterised by mg = as in the 
anomalous phase and exhibit a length scale that grows 
with the time that has elapsed since a quench from above 
Tc- It can be readily shown that if multiple solutions 
to the equations of motion exist with different activity, 
then the field s acts to select the solutions with the larger 
(s > 0) or smaller (s < 0) activity, leading to a transition 
at s = 0. However, while the aging dynamics of the p ~ 2 
spin-glass of ([T]) can be solved [251 by a similar method 
to that of Cugliandolo and Dean [20|, we have not yet 
established any clear connection between these dynam- 
ics and the anomalous disordered phases discussed here. 
This too remains an area for future study. 

We also compare the results shown here with those 
obtained for kinctically constrained models [l^ . In both 
cases, active and inactive phases coexist at the s = axis. 
However, there are two important differences. Firstly, 
in the kinetically constrained models considered in p^ . 
the transition is first-order, signalling the coexistence of 
active and inactive solutions to the equations of motion. 
These have been interpreted as "ergodic fluid" and "non- 
crgodic glass" states 0, H, [H, [11] . On the other hand, 
the continuous transition in the soft spin models is second 
order: the anomalous phase is characterised by a diverg- 
ing correlation length that we have tentatively attributed 
to the growing length scale associated with the aging be- 
haviour of the system. Taking the large-time limit of the 
aging solution, the activity of the system approaches that 
of the ordered state: the active (aging) and inactive (or- 
dered) phase are not separated by a gap in the activity, 
unlike the kinctically constrained models. 

Secondly, wc emphasise that in the kinetically con- 
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strained models, the s = axis of the phase diagram 
belongs to the active phase. Introducing any s > leads 
immediately to an inactive phase that is qualitatively dif- 
ferent from the unbiased state at s = 0. On the other 
hand, in the models considered here, the s ~ axis be- 
longs to the inactive (ordered phase): it is the introduc- 
tion of any s < that leads to an active phase that differs 
from the unbiased steady state. 

Finally wc note that contrary to its p > 3 coun- 
terparts, the thermodynamic properties of the p = 2 
spin glass have only a single transition temperature and 
do not display any kind of replica symmetry breaking. 
Here, we have analysed this problem by diagonalising 
the quadratic dynamical action: a method that applies 
only for p = 2. However, the same results can be verified 
using the replica trick and integrating out the disorder. 
The application of such methods to models with p > 3 
would provide further insight into the behaviour of large 
deviations of the activity in 'glassy' models. 
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So{t) = / dt' [X{t')X{t') + X{t')X{t') + AuX{t')X{t')] 
Jo 

(A4) 



and 



Jo 



— + X{t') 13 J, 



+ m <t>i~<i>i\ (A5) 



where we again omit the dependence of the fields (f) and 
(j) on the time t' , for brevity. 

In the N oo limit, the integrals over X, X, A and A 
can be carried out through a saddle point approximation. 
We replace X{t), X{t), X{t) and X{t) by their saddle-point 
values X, Xi ^* f-iid A* . In particular, differentiating the 
action (iV^o -1-5*1) with respect to X and X, we arrive at 
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A* = -4ux 
X* = -Aux . 



(A6) 
(A7) 



Thus, performing the saddle-point integrals in (jA3[) leads 
to dni) in the main text. Similarly, differentiating with 
respect to A and A leads to the self-consistency equations 
(HH) and (EOl). 



APPENDIX A: SADDLE-POINT INTEGRATION 
OF THE DYNAMICAL ACTION 

Here, we show how the functional integral can be 
cast in the form (|17p . We define 



(Al) 



which is proportional to the growth rate of the activity 
K defined in pT|) . and the related quantity: 



(A2) 



Both X and X are self-averaging quantities in the ther- 
modynamic limit of large N. It is therefore convenient 
to constrain these quantities with Lagrange multipliers 
and then to integrate over these constrained quantities 
by a saddle-point method. Writing the (time-dependent) 
Lagrange multipliers as A, A, Eq. (|15p becomes: 



Z{s,t) 



vxvx^^vm 

2111 27ri 



-WSo(t)-Si(t) 



(A3) 



APPENDIX B: COMPUTATION OF x IN THE 
LOW-TEMPERATURE DISORDERED PHASE 



Here we discuss the solutions of the self-consistency 
equations ^ and pi)) f or s < and T < T^. Using 
the notation of Sec. IIIIEl we can write the integrals of 
Equs. (gT]) and (gS]) as 



h = 



h = 



4+1- 



dz ^4(x' ~ z) ~ (x' - z)2 
' 27r + y' 

dz -2^/4(2;' - z) - {x' - zY 



, 2tt 



^Jz^ 



v 



(Bl) 
(B2) 



Since wc are at (3 > (3c and s < 0, we have x',y' > 0. 
The small s limit becomes the limit ?/' —s- but if we 
take the limit for y' — > keeping a;' > we have 
Ii{x',y' ^ 0) ^ 00, and the self-consistency condition 
cannot be satisfied. We therefore take both x' and y' to 
zero together: we assume that y' <C x'^ which can be ver- 
ified a posteriori through the solution to (jSip . The result 
is x' ~ [ln(l/y')]~^, consistent with our assumptions. 
We start by splitting the integral in (jBip into three 
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parts: 



^ dz v^x'(4-a;') 
-X' 27r y^y' + z2 



^4(a;' - z) - {x' - zY - ^x'(4-x') 

yz2~r7 

dz ^\{x' - z) - (.x' - z)2 



27r 



(B3) 



In the limit of y' <C x'^, the first integral has a diver- 
gent contribution {\fx! jii) In [Ax'"^ /y'). In the second and 
third parts we can take the limit y' ^ directly in the 
integrand. If we then take the limit of small x' then we 
find that the second integral vanishes as 0{\/x') while 
the third approaches unity. Thus we arrive at (|49)) . 

We now evaluate the {y' /x' ) — > limit of the expres- 
sion (jB2p . We use a similar method, spliting the integral 
into two terms: 



In the first integral we introduce w = zjx' and a 
y' jx' , thus arriving at 



(B5) 

The leading behaviour of this quantity can be evaluated 
by setting directly a = x' = Q into the integral, so that 
this contribution to h is -4(x')3/2/(37r)[l + o(l)]. 



For the second integral in Eq. (jB4[) . we have z^ > x''^ 
so we expand the integrand in powers of cr = y' I z^ . The 
leading term vanishes and the second term is 0{y' /^fx'). 
Writing {y' / ^/x) = {x')^^'^{y' /x'^), this term is smaller 
first term in (|B4p which is 0{{x')^^^). Thus we arrive at 



h = 1- 



' ^^4{x'-z)-{x'~zyl^=L= 



— ^ rj-, 

+ / —^J^{x'~z)-{x'-zf 



27r 



yjz'^ + y' 



li « 1 



4(x') 



3/2 



(B^) as given in Equ. ([501) of the main text. 



(B6) 
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